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Abstract: We suggest a way of improving the probes on dimension -6 CP-conserving HVV 
interactions {V = W, Z, 7 ), from the LHC data on the Higgs boson to be available in the 
14 TeV run with an integrated luminosity of 3000 fb“^. We find that the ratios of total 
rates in different channels can be quite useful in this respect. This includes ratios of event 
rates in (a) different final states for the Higgs produced by the same production mechanism, 
and (b) the same final state from two different production modes. While most theoretical 
uncertainties cancel in the former, the latter helps in the case of those operators which 
shift the numerator and denominator in opposite directions. Our analysis, incorporating 
theoretical, systematic and statistical uncertain, leads to projected limits that are better 
than the strongest ones obtained so far from precision electroweak as well as LHC Higgs 
data. Moreover, values of the coefficients of the dimension -6 operators, which are allowed 
in disjoint intervals, can have their ranges narrowed down substantially in our approach. 
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1 Introduction 

The ATLAS and CMS experiments at the Large Hadron Collider (LHC) have discovered 
a neutral spinless particle that closely matches the description of the Higgs boson [1, 2] 
which is responsible for masses of elementary particles, according to the standard model 
(SM) of electroweak interactions. While this ties the final knot on the framework embodied 
in the SM, there are many reasons to believe that there is more fundamental physics at 
higher energies. The reason for such expectation can be traced to many issues, including 
the unexplained replication of fermion families, the source of dark matter in the universe, 
and the problems of naturalness and vacuum stability involving the Higgs boson itself. The 
Large Hadron Collider (LHC) has not revealed any direct signature of new physics so far. 
However, one is led to suspect that such physics should affect the interaction Lagrangian 
of the Higgs boson. This generates, for example, effective operators of dimension -6 con¬ 
tributing to HVV interactions, with V = W, Z, 7 . Probing such effective couplings for the 
recently discovered scalar is therefore tantamount to opening a gateway to fundamental 
physics just beyond our present reach. 

Such ‘effective’ interaction terms better be SU{2) x 77(1) invariant if they arise from 
physics above the electroweak scale. Constraints on such terms have already been studied, 
using precision electroweak data as well as global fits of the current Higgs data [3-37]. 
Recently, CMS has published an exhaustive study on anomalous HVV couplings [38]. 
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Many studies have considered anomalous Higgs couplings in context of future lepton col¬ 
liders [39-44], The general conclusion, based on analyses of the 8 TeV data, is that several 
(though not all) of the gauge invariant, dimension-6 HVV terms have been quite strongly 
constrained by the EW precision and LHC data (as discussed in section 3) [3-37], It still 
remains to be seen whether such small coefficients can be discerned with some ingeniously 
constructed kinematic distributions. Some work has nonetheless been done to study such 
distributions [45-49], in terms of either the gauge invariant operators themselves or the 
structures finally ensuing from them. At the same time, it is of interest to see if meaning¬ 
ful constraints do arise from the study of total rates at the LHC. The essence of any probe 
of these anomalous couplings, however, lies in pinning them down to much smaller values 
using the 14 TeV runs, as common sense suggests the manifestation, if any, of new physics 
through Higher Dimensional Operators (HDO’s) with small coefficients only. 

We show here that the relative rates of events of different kinds in the Higgs data 
can allow us to probe such effective interactions to levels of smallness not deemed testable 
otherwise [50, 51]. This happens through (a) the cancellation of theoretical uncertainties, 
and (b) the fact that some ratios have the numerators and denominators shifting in opposite 
directions, driven by the additional interactions. Thus the cherished scheme of finding 
traces of new physics in Higgs phenomenology can be buttressed with one more brick. 

We organise our paper as follows: we summarise the relevant gauge invariant operators 
and the interaction terms in Sec. 2. In Sec. 3, we introduce three ratios of cross-sections 
as our observables. The results of our analysis are explained in Sec. 4. We summarise and 
conclude in Sec. 5. 

2 Higher dimensional operators 

In order to see any possible deviations from the SM in the Higgs sector, we will follow the 
effective field theory (EFT) framework. We consider SU{2)lxU (l)y invariant operators of 
dimension up to 6, which affect Higgs couplings to itself and/or a pair of electroweak vector 
bosons. While a full list of such operators are found in [52-55], we have concentrated here 
on dimension-6 CP-conserving operators which affect Higgs phenomenology. They include: 

• Operators which contain the Higgs doublet and its derivatives: 

0$,i = (D^$)t$$t(D'^^); 0$,2 = ^dM($^d>)d'‘($^^); = ^(^^^)" (2.1) 

• Those containing $ (or its derivatives) and the bosonic field strengths : 

Ogg = Obw = Oww = 

= = OB = (D^4>)tH/^-(D,$), (2.2) 

where 

= i ^aaW^ = i - 

2 2 
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and g, g' are respectively the SU(2 )l and U(1 )y gauge couplings. — dyW^ — 

g^abcwb^^c^ = d^By - dyB^ and = d^Gl - dyG^ - gs^^^G^Gl. The covariant 
derivative of is given as = [dfj, + ^g'B^+ig^W^)^. The Lagrangian in the presence 
of the above operators can be generally expressed as: 

Cz^ + (2.3) 

^ ' i 

where in addition to the dimension -6 (D 6 ) operators, we also allow for the SM-like HWW 
and HZZ couplings to be scaled by a factor k. While k 7 ^ 1 is indicative of certain kinds 
of new physics, we are specially interested in this study in the new observable features 
associated with the HDOs. Therefore, we have set k = 1 for simplicity.^ 

No operator of the form Ogg is assumed to exist since we are presently concerned 
with Higgs interactions with a pair of electroweak vector bosons only. The operator 
is severely constrained by the T-parameter (or equivalently the p parameter), as it alters 
the HZZ and HWW couplings by unequal multiplicative factors. As far as HZZ and 
HWW interactions are concerned, only scales the standard model-like couplings (k), 
without bringing in any new Lorentz structure. This amounts to a renormalization of the 
Higgs field. It also alters the Higgs self-coupling, something that is the sole consequence 
of as well. 

In view of the above, we focus on the four operators Oww, Gbb, Gw and Ob- We 
do not include the operator Obw = Bf^yW^'^^ in the present analysis, because it mixes 

the Z and 7 fields at the tree level, violates custodial symmetry (by contributing only to 
the Z-boson mass) and is, therefore, highly constrained by the S and T-parameters at the 
tree level [4]. The effective interactions that finally emerge and affect the Higgs sector are 

^eff = gilw {W+yW-^d'^H + h.c.) + g%^^^ HW+yW-^^'^ 

+ < 7 $;^ Z^yZi^d^'H + g%\z HZ^.yZ^'^ 

+ A^yZ^^d’^H + HA^yZ>^- + gn^-.HA^yA^'^, (2.4) 

where 

f 9Mw\ fw ( 2 ) _ f gMw\ „ 

1 ^2 ) 2 ’ ~ i a 2 y 

(gMw\ c^fw + (2) _ (gAIw\ s^/bb + c^fww 

V A2 2c2 ’ ^HZZ- \ ) 2c2 

(gAlw\ s{fw - Ib) ( 2 ) _ f gMw\ s{s^fBB - <^fww) 

\ ) 2 c ' ’ “ V A" y c 

/ gAIw 
V A 2 

with s (c) being the sine (cosine) of the Weinberg angle. Besides, the operators Ow-, Ob 
and Owww also contribute to the anomalous triple gauge boson interactions which can be 

^Possible constraints on the departure of k from unity have been obtained in the literatnre from global 
fits of the Higgs data (See for example [3-36]) 


s^Ubb + fww) 


(2.5) 


9hww ~ 


9hzz = 


^HZ'y 


dH'j'y 


- 3 - 















summarised as 


Hwwv = -igwwv l^r {w;,W-^V- - + KyW+W-V^^^ + 

( 2 . 6 ) 


where gwW'y = g s, gwwz = 5 c, Ky = 1 + Any and gf = I + Agf with 


/lf2 

ifw + Ib) : 


o _2 ^2 


\ 


(2.7) 


The already existing limits on the various operators discussed above are found in 
numerous references [3-6, 10]. Even within their current limits, some of the operators are 
found to modify the efficiencies of the various kinetic cuts [9, 14]. The question we address 
in the rest of the paper is : can these limits be improved in the next run(s) through careful 
measurement of the ratios of total rates in different channels? As we shall see below, the 
answer is in the affirmative. 


3 Ratios of cross-sections as chosen observables 

The four HDOs under consideration affect Higgs production as well as its decays, albeit 
to various degrees. For example, HDO-dependent single Higgs production processes are in 
association with vector bosons {VH) i.e. pp ^ VH (where V = {W, Z}) and vector-boson 
fusion (VBF). We show the production cross-sections in these channels at 14 TeV in Fig. 1, 
as functions of the four operator coefficients (/*) taken one at a time.^ The relevant decay 
channels which are dependent on such operators are H —)• WW*, ZZ*,'y'y, Zj. Fig. 2 con¬ 
tains these branching ratios (HR) as functions of the four coefficients under consideration. 

The VBF and VH rates are sensitive to fww and fw, but depend very weakly on 
Jbb and fs, while the cross-section a{pp —>• WH), is completely independent of /bb and 
/b)- The HDO effects in 77 —)• 77 and H —)• Z 7 for fi ~ 0(1} ^ is of the same order as 
the loop-induced SM contribution unlike in the case of the HWW and HZZ couplings. 
Therefore, becomes highly sensitive to fww and Jbb- Consequently, the 7-1-8 

TeV data already restrict their magnitudes. Bounds on all these operators in a similar 
framework can be seen in Table VI of Ref. [4] and also in Ref. [3] . In Ref. [4] , the bounds 
have been presented at 90% CL by varying multiple operators at the same time. These 
bounds have been obtained by considering the LHC data as well as constraints from on 
the oblique parameters, viz., S,T and U. Bounds coming from the oblique parameters are 
generally weaker than those obtained from the LHC data as can be seen in Ref. [3]. These 

^We have used CTEQ6L1 parton distribution functions (PDFs) by setting the factorization (fip) and 
renormalization scales {v-r) at the Higgs mass {Mr = 125 GeV). 

®If the operators arise from loop-induced diagrams which imply ‘loop factors’ in denominators of the 
effective interactions, 0(1) TeV“^ coefficients imply strongly coupled theories [11, 56]. However, if such 
operators originate from tree-level diagrams, then 0(1) TeV“^ coefficients imply weakly-coupled theories. 
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(c) 


Figure 1. Higgs production cross-sections for the VBF and VH channels in presence of HDOs at 
14 TeV. Here the operators are varied one at a time. 


limits may not be applicable when the analysis is performed varying one operator at a 
time. 

Based on the above information, we set out to hnd observables which are sensitive to 
/j < 5 TeV“^ in the High luminosity run at the LHC. It is not completely clear yet how 
much of statistics is required to probe such small values with various event shape variables. 
On the other hand, the more straightforward observables, namely, total rates in various 
channels, are also fraught with statistical, systematic and theoretical uncertainties which 
must be reduced as far as possible when precision is at a premium. 

An approach that is helpful is looking at ratios of cross-sections in different channels. 
In this paper, we invoke two kinds of ratios. First, we take ratios of events in two different 
final states arising from a Higgs produced via the same channel (in our case, gluon fusion). 
Such a ratio enables one to get rid of correlated theoretical uncertainties (CThU) such as 
those in PDF and renormalisation/factorisation scales. They also cancel the uncertainty in 
total width which is correlated in the calculation of BRs into the two hnal states. Secondly, 
we consider the ratio of rates for the same final state for two different production channels 
(such as VBF and VH). Although the uncertainty in the BR cancels here, the theoretical 
uncertainties at the production level do not. Moreover, since the final state is same in this 
case, some systematic uncertainties which are correlated (related to identification, isolation, 
trigger etc.) will also get cancelled. However, this is helpful in another manner. For some 
of the operators, the /j-dependent shifts with respect to the SM are in opposite direction 
for the numerator and the denominator in such ratios. The result is that the net deviation 
adds up, as shown in subsection 3.2. We shall see that the use of both these kinds of 
ratios (including those involving the channel can capture the HDO coefficients at a 
level unprecedented, going down to values where new physics can show up. 

3.1 Observable sensitive to Oww Obb' Hi 

As has been noted earlier, BRj 7 _>..y.y (Fig. 2c) is highly sensitive to two of the operators, 
namely, Obb and Oww- Therefore, we propose to probe them in the 77 channel, with the 
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Figure 2. Branching ratios oi H ^ WW*, ZZ *in presence of HDOs. The operators are 
varied one at a time. 


Higgs produced through gluon-gluon fusion (ggF). This final state is clean for reconstruc¬ 
tion, and has high statistics. We should mention here that if we consider the simultaneous 
presence of more than one operators, then there is a “blind-direction” in the parameter 
space fww ~ —fsB where mimics the SM value. This is because the higher¬ 

dimensional part of the vertex is proportional to fww + Jbb- Also, for the non-trivial 
range fww = fBB ^ —3, BRh^'y'y mimics the SM value, due to parabolic dependence of 
the diphoton rate on the HDO coefficients. Therefore, the Higgs produced through ggF 
followed by its decay to 77 cannot be used alone to probe these two ‘special’ regions of the 
parameter space. We construct the observable 


niifi) 


(TggF X 

CJggF X BRH ^ WW *^ 2 e 2 u { fi ) 


(3.1) 
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where i = e, fj, and /j’s are the operator coefficients. As explained earlier, the CThU in 
production as well as total width cancels here; so does the AT-factor in the production rate. 
Clearly, TZi can also be expressed as the ratio of two signal strengths as follows, 


AffU) ^ (UggF X BRn^-y-yf^ 

{(^ggF x'BRH^WW*^2£2uf^ 

Therefore, already measured 77 and WW* signal strengths can be used to constrain 
the operator coefficients affecting the ratio TZi. The efficiency of acceptance cuts does 
not affect the results, for values of fww and fsB which are of relevance here because for 
such small values of the parameter coefficients the change in experimental cut-efficiencies 
is negligible. On top of that, for the ggF production mode, these operators only affect 
the decay vertices and hence the cut-efficiencies are but modified by a very small extent. 
We must also note that in defining TZi a full jet-veto (0-jet category) has been demanded 
for both the numerator and the denominator to reduce the uncertainties related to the 
different jet-requirement in the hnal state. Besides, in the denominator, the WW* pair is 
considered to decay into both same flavour (ee -|- fi/j.) and different flavour {eg. + ge) final 
states to improve the statistics. 



3.2 Observable sensitive to Oww and Ow- ^2 

It turns out that the fww and fw affect (to one’s advantage) the ratio of events in a 
particular Higgs decay mode in the VBF and VH channels. This captures the new physics 
at the production level. By considering the same final states from Higgs decay, some 
theoretical uncertainties in the decay part cancels out. The production level uncertainties, 
including the AT-factors, however, do not cancel here. In our calculation, the next-to-next-to 
leading order (NNLO) /C-factors have been assumed to be the same as in the SM, expecting 
that the presence of HDO does not effect the AT-factors much. For precise estimate of the 
observed ratio, one of course has to incorporate the modified cut efficiencies due to the new 
operators, though such modifications may be small. The other, important advantage in 
taking the above kind of ratio is that, for not-too-large fww or fw (in the range [—5, -|-5]), 
the deviations of the VBF and VH cross-sections are in opposite directions. The generic 
deviation for the rate in any channel can be parametrized as 

HDO \y {‘i _i_ ^ ^ fQ q\ 

‘^prod. ~ '^prod. X (1 + Oprod.j • (o-o) 

From Fig. la, dvBF is positive in the range fww, fw > 0. On the other hand, in the same 
region of the parameter space, 5 vh is negative as evident from Figs, lb and Ic. Hence, on 
taking the ratio /o'y^^, the deviation from SM is 

X (l + <5 vbf — 5vh + C>((5^)) . (3.4) 

O-VH cr^H 

Thus this ratio further accentuates the deviation from SM behaviour. As an example, if 
we consider the parameter choice fw = 2, then <5 vbf ~ 3.6% and (5 wh ~ 10%. However, 
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from the ratio, the combined (5vbf+wh ~ 15%, which is a clear indication of why we should 
consider such ratios. We thus define our next observable 

^ o'vbf(/0 X . . 

o'wnifi) X X BRw-s.fl/’ 

where the 77 final state has been chosen because of its clean character and reconstructibility 
of the Higgs mass. It should be remembered, however, that fww,fBB in the range —3 to 
0 causes the diphoton branching ratio to undergo a further dip. This can adversely affect 
the statistics, and thus the high luminosity run is required for an exhaustive scan of the 
admissible ranges of the above coefficients. 


3.3 Observable sensitive to Ob'- ^3 


The operator Ob is sensitive to H ^ ZZ* and H —)• Z 7 . In the former mode, the 
sensitivity of /b is limited (see the green curve in Fig. 2b) and can be appreciable only for 
larger Jb- The partial decay width Th^Z')) on the other hand is rather sensitive to all the 
four operators under study (Fig. 2d), primarily due to the fact that the new FfZy vertex 
contributes practically as the same order as in the SM. However, the present statistics in 
this channel is poor [57, 58]. We expect better bounds on Oww-, Obb and Oy/ from the 
measurements of and 7 ^ 2 - We use TZ^ for the 14 TeV 3000 fb“^ run to constrain Jb 
only, for which other channels fail. In the same spirit as for TZi, we thus define our third 
observable 


n^Ui) = 


CJggF X 


(3.6) 


CJggF X BRj^_^wTV*-5>2f2y(/i) ’ 
where i = e, n and here again the CThU cancels. Here also, we must note that in defining 
7^3 a full jet-veto has been demanded for both the numerator and the denominator. For 
the numerator, the Z boson’s decay to both an electron pair and a muon pair is considered. 
Besides, in the denominator, the WW* pair is taken to decay similar to the TZi case. 

Comparison with the K-framework : In principle, studies in terms of ratios in 
different channels can be carried also in the K-framework [ 8 , 59-61] in which couplings are 
modified just by scale factors. It should, however, be remembered that the present analysis 
involves new Lorentz structures and hence brings non-trivial interference terms in the 
squared amplitudes. Unlike the situation with overall scaling, this prevents the cancellation 
of the modifying couplings when one considers ratios of events taking (SM-I-BSM) effects 
into account. 

Even though the ratio 7^i (T^s), dominated by 7/77 (HZ^) vertex, contains no new 
Lorentz structures, it is still sensitive to the HDOs due to the presence of the HWW 
vertex in the denominator. Therefore, these ratios, although apparently similar to ratios 
employing the k- framework, are different in practice. 7^2 is a ratio of (Jvbf and ctwh which 
are sensitive to the operator coefficients as shown in Fig. 1. In the k- framework, avBF is 
dominated by the WWH vertex and hence nww will approximately cancel in 77 - 2 . On the 
other hand, there will be no trivial cancellations between the numerator and denominator 
in the HDO-framework. 
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4 Results of the analysis 


For our subsequent collider analysis, the chain we have used is as follows - first we have 
implemented the relevant dimension-6 interaction terms as shown in Eq. (2.4) in Feyn- 
Rules [62], and generated the Universal FeynRules Output (UFO) [63] model files. These 
UFO model files have been used in the Monte-Carlo (MC) event generator Mad- 
Graph [64] to generate event samples. Next, the parton-showering and hadronisation are 
performed using Pythia [65] and finally detector level analyses is carried using Delphes [66] 
Before we discuss the phenomenological aspects of the aforementioned observables, we 
re-iterate below the various kinds of uncertainties considered. The two major classes of 
observables where these uncertainties arise are as follows: 

• Same production channel but different final states: 

In such cases (as in TZi and T^s), the correlated uncertainties lie in PDF-l-as, QCD- 
scale and in the total Higgs decay width, Fjif. However, uncertainties in the partial 
decay widths are uncorrelated Statistical uncertainties for distinct final states are 
always uncorrelated and are retained in our analysis. We also assume some systematic 
uncertainties, whenever shown, to be fully uncorrelated. All surviving uncertainties 
are added in quadrature to estimate total uncertainties related to our observables. 

• Different production channels but same final state: 

For such observables (7^2 in our definition), the only correlated uncertainty is in 
All other uncertainties are uncorrelated and hence are added in quadra¬ 
ture (including the uncertainties in the numerator and the denominator of the ratio 
7 ^ 2 )- Beside the already mentioned theoretical uncertainties, we also encounter some 
additional theoretical uncertainty related to the QCD-scale in the WH mode, which 
we separately discuss in subsection 4.3. 

We further assume that the percentage uncertainties remain same even after the inclu¬ 
sion of the anomalous couplings. In order to illustrate, how the uncertainties are taken into 
consideration, we list the theoretical uncertainties related to relevant Higgs BR and total 
width in Table I, and related to various production cross-sections in Table 2. In Table 3, 
we present the number of surviving events after the selection cuts in the SM at 14 TeV 
with 3000 fb“^ integrated luminosity in the pure production modes. These numbers are 
taken from Refs. [ 68 , 69] except for the 77 channel in the VBF production mode, which 
we estimate by applying a fixed pT-cut (keeping other cuts are same as in Ref. [ 68 ] ) of 50 
GeV on both the tagged jets instead of r/-dependent jet selection cuts as used in the same 
reference. The number of events have been computed by removing the contaminations from 
other production mechanisms which will reduce the number of events and hence enhance 
the statistical uncertainties (which roughly goes as \/Ns + Nb/Ns, with Ns and Nb being 
respectively the number of surviving signal and background events after selection cuts). 

^We must mention here that and Th^zj have tiny correlations with Th^ww because of the 

IP-boson loop in the former two cases. However, in this present analysis we neglect such small correlations 
and consider these partial decay widths to be mostly uncorrelated 
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SM Quantity 

Value 

+ve uncert. % 

—ve uncert. % 


2.28 X 10-3 

-b4.99 

-4.89 

BRh^ww* 

2.15 X 10"^ 

-b4.26 

-4.20 


1.07 X 10-^ 

-b0.16 

-0.16 


1.06 X 10-^ 

-b0.15 

-0.15 

BRi7-s>Z7 

1.54 X 10-3 

-b9.01 

-8.83 


3.36 X 10-2 

-b0.004 

-0.004 


3.37 X 10-2 

-b0.007 

-0.007 

Total Tjy 

4.07 MeV 

-b3.97 

-3.94 


Table 1. i BRz->.« and total Higgs width F/f (MeV) 

and their % uncertainties (+ne and —ve refer to positive and negative uncertainties respectively) 
for a Higgs of mass 125 GeV {mw = 80.385 GeV and mz = 91.1876 GeV). These 
numbers are taken from the LHG Higgs Cross Section Working Group page [67]. 


Process 

a (pb) 

-hQCD-Scale % 

-QCD-Scale % 

+(PDF+as) % 

-(PDF+as) % 

ggF 

49.47 

+7.5 

- 8.0 

+7.2 

- 6.0 

VBF 

4.233 

+0.4 

-0.5 

+3.3 

-3.3 

WH 

1.522 

+ 0.8 

- 1.6 

+3.2 

-3.2 

ZH 

0.969 

+4.0 

-3.9 

+3.5 

-3.5 


Table 2. The cross-sections of relevant Higgs production {mn = 125 GeV) channels and their 
QCD-Scale and PDF-t-Us uncertainties in %. These numbers are again taken from the LHC Higgs 
Cross Section Working Group page [67]. 




7^2 

^3 

jS^num 

47724 (77 in ggF) 

194 (77 in VBF) 

1989 (Z 7 in ggF) 

Ajnum 

3.16 X 10® 

1041 

691931 

Nden 

40850 {WW* in ggF) 

238 (77 in WH) 

40850 {WW* in ggF) 

Nden 

366450 

995 

366450 


Table 3. Number of surviving events (taken from Refs. [ 68 , 69]) after the selection cuts in the 
SM at 14 TeV with 3000 fb“^ integrated luminosity. These numbers are used to compute the 
statistical uncertainties (which goes as \/Ns + Nb/Ns, where Ns and Nb are respectively the 
number of surviving signal and background events after selection cuts) related to the numerator 
and denominator of the three observables. Number of events in the VBF ( 77 ) channel is computed 
by applying a hxed pr-cut (keeping other cuts same as in Ref. [ 68 ] ) of 50 GeV on both the tagged 
jets instead of ry-dependent jet selection cuts as used in the same reference. Number of events for 
77 in TZi, in 7?.3 and WW* for 'R.^ and 'R.-i are obtained after putting 0-jet veto and demanding 
only ggF events. The superscripts num and den signifies the numerators and denominators of the 
three observables. 

For instance, the reported number of 77 events for an integrated luminosity of 3000 fb“^ 
is 49200 with a 3% contamination from VBF (Table 3 in Ref. [ 68 ]). In our analysis we 
have used Ns = 47724 (= 0.97 x 49200) to compute the statistical uncertainty. Similarly 
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TZi 

TZ2 

^3 

2.87 % 

13.83 % 

29.63 % 


Table 4. Statistical uncertainty for the observables T^i, 0-2 and The numbers are obtained 
after doubling the number of signal and background events given in Table 3 in order to account for 
both ATLAS and CMS experiments. 



7^l 

7^2 

TZ3 

Numerator 

2.5% (77 in ggF) 

9.1% (77 in VBF) 

3.1% (Z 7 in ggF) 

Denominator 

3.4% {WW* in ggF) 

5.0% (77 in WH) 

2.8% (WW* in ggF) 


Table 5. Systematic uncertainties used in our analysis to compute the total uncertainties related 
to the three observables. The numbers shown here are combination of various types of relevant 
systematic uncertainties added in quadrature taken from Refs. [57, 70, 71]. 

a 30% contamination in the VBF category due to ggF (Table 3 in Ref. [ 68 ]) has also 
been taken into consideration. In doing so, we are giving conservative estimates on the 
statistical uncertainties. All entries in Table 3 are shown after removing contamination 
to compute conservative statistical uncertainties. We must note that, while computing 
the statistical uncertainties (as shown in Table 4) for all the three ratios, we double the 
number of events in Table 3 to roughly accommodate two independent experiments to be 
performed by ATLAS and CMS. Here, we assume that ATLAS and CMS will analyse the 
same channels with similar set of selection cuts and will roughly obtain same number of 
events in the actual experiment. It is also assumed that the overall performance of ATLAS 
and CMS will be similar, integrated over a large luminosity. In future, when the data be¬ 
come actually available, one would be able to compute the exact statistical uncertainties. 
However, we must note that one should actually take the number of events in the side-band 
{^side-band) in Order to compute the statistical uncertainties. The procedure we follow 
gives conservative values for the statistical uncertainties. In future, the actual experiments 
will provide us Nside-band which will allow us to compute accurate statistical uncertainties. 
However, the side-band analysis is beyond the scope of this paper as the data for the 14 
TeV run at 3000 fb“^ is yet unavailable. 

We also use some systematic uncertainties in our analysis as listed in Table 5 (Refs. [57, 
70, 71]). In the future, it is quite expected, various systematic uncertainties will reduce by 
improving their modelling. To be conservative, we have used various important uncorre¬ 
lated systematic uncertainties as used in Refs. [57, 70, 71] for 7-1-8 TeV analysis. For the 
observable TZi, since we are applying same jet veto {i.e. 0-jet category), the systematic 
uncertainties related to the jet energy scale, jet vertex fraction etc. will not be present. 
On the other hand, due to the different final state, systematic uncertainties related to the 
photon and lepton identification and isolation, missing energy trigger etc. will remain. In 
a similar fashion, for 7^2 and TZ^ various correlated systematic uncertainties will cancel 
between their respective numerator and denominator. 

Next, we consider the ratio TZi in the light of both the existing data and those predicted 
for the high energy run. For 7^2 and TZ^, only a discussion in terms of 14 TeV rates is 
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relevant, as the currently available results have insufficient statistics on these. 


4.1 7^l @ 7+8 TeV 

Before predicting the bounds from the 14 TeV HL run, let us form an idea about the 
constraints from the 7+8 TeV Higgs data in the 77 and WW* channels. In Table 6 , 
we show the exclusive signal strengths in the 77 and WW* hnal states through the ggF 
production mode as reported by ATLAS [70, 71] and CMS [72, 73]. 

We must emphasize that the categorization introduced by the ATLAS and CMS ex¬ 
periments are used to enhance the sensitivity for the Higgs boson signal (Tables H and 
HI in Ref. [70]). The signal strengths {g) shown in Fig. 17 include these contaminations. 
These signal strengths are further combined to give specific production categories as shown 
in Fig. 18. For instance g for ggF categories is the combination of the four categories, viz. 
central low Pt^, central high Pt^, forward low Pxt and forward high Pt^. Therefore, the 
g for specific categories in Fig. 18 is not exclusive. However, while obtaining the g for 
a specific production mode in Fig. 19, the effect of contaminations are properly removed 
(by knowing the amount of contaminations from Monte-Carlo simulation for the SM) and 
therefore, these are the exclusive signal strengths. The removal of contaminations includes 
not only the subtraction of production mechanisms that are not of interest but also the 
propagation of errors. The experiments have taken into account the impact on the statis¬ 
tical, systematic and theoretical errors for the extraction of the exclusive signal strengths. 
Therefore, the exclusive g will generally contain larger uncertainty. For example one can 
see that the error on the global signal strength is signihcantly better than that extracted 
for individual production mechanisms. For instance, in Ref. [70], where ATLAS reports on 
signal strengths with the di-photon channel, the global signal strength is + = 1.17 ± 0.27, 
which leads to an accuracy of 23%, whereas for the signal strength of gluon-gluon fusion 
(ggf) gggf = 1.32 + 0.38, corresponding to an accuracy of 29%. Same applies to the results 
reported by CMS in Ref. [72]. 

Here we statistically combine the signal strengths for a particular final state as reported 
by the two experiments, using the following relations 


— = y"—■ 

^ af ’ 

7, ^ 


g _ gi 
+2 ^ aj ’ 

7. 


(4.1) 


where a {g) refers to the combined Icr uncertainty (signal strength) and cjj {gi) signifies 
the corresponding uncertainties (signal strengths) in different experiments. 

We compute all the surviving correlated theory errors and subtract them in quadrature 
from the errors in the numerator and denominator of the ratio TZi, viz. = g^^^^ x 

((TggF X and x (uggp x BR^_,.rvty*)^^ x y~]^ 

Fig. 3, the red line is the theoretically computed which is independent of the centre of 
mass energy since TZi is actually a ratio of two BRs. The outer (light green) band shows 

®For instance, the error associated with combined (ATLAS-I-CMS) —>■ 77 ) i.e. ±0.26 consists of 

theoretical, statistical and systematic uncertainties and, by subtracting the CThU (±0.13) in quadrature 
we get (±0.22) which will finally contribute to the uncertainty related to the numerator of TZi. 
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Experiment 

fj,{H 77 ) in ggF 

WW* 21$b) in ggF 

ATLAS (@ 7+8 TeV) 

1 00+0.38 
-‘-•'^^-0.38 

-‘-•n^-0.26 

CMS (@ 7+8 TeV) 

1 19+0.37 
-‘-•-‘■^-0.32 

n- ‘a o.23 

Combined 

1.21 + 0.26 

0.88 + 0.19 


Table 6 . Measured Higgs Signal strengths in the 77 and WW* modes where Higgs is produced 
through only ggF channel using -/s = 7 + 8 TeV data by ATLAS [70, 71] and CMS [60, 72]. Here 
we have combined the ATLAS and CMS signal strengths for a particular final state and production 
mode using Eq. 4.1. 


the uncertainty comprising of the uncorrelated theoretical, statistical and systematic parts 
and the inner (dark green) band represents the total uncorrelated theory uncertainty. The 
black dashed line gives the experimental central value of TZi. The ratio, 7Zi is almost 
completely dominated by (since 'SRh^ww* is not so sensitive on HDOs) and 

therefore highly sensitive to the operators Oww and Obb- The parabolic nature of the 
BRn^-y'-/ as functions of fww and /bb leads to two disjoint allowed ranges of fww = /bb ~ 
[—3.32, —2.91] U [0.12,0.57] as shown in Fig. 3. We should mention that the region between 
these two allowed ranges shows extremely low values of because of destructive 

interference between the SM and HDO might leads to poor statistics. If both Oww and 
Obb are present simultaneously with almost equal magnitude and opposite signs, the 
observable TZi closely mimics the SM expectation, and to probe that ‘special’ region of 
parameter space we need to go for other observable like 7^2- The operators Ow and Ob 
are mostly insensitive to this observable mainly because BR.^.-,, is independent of these 
operators and the dependence of BRu/pp* on all four operators is comparatively weak (see 
Fig. 2a) 

We compare our results with the existing bounds on these operators as obtained in lit¬ 
erature. For instance, the limits obtained in Fig. 3 (left panel) of Ref. [3] on Oww and Obb 
at 68% CL are [—3.23, —2.61] U [—0.35,0.27] (in TeV“^) for the ATLAS case. In obtaining 
these limits, they varied one operator at a time. This is similar in approach to our study 
where we have given a framework where one operator is varied at a time. Our bounds are 
in very good agreement with their results. The slightly different limits obtained by us are 
due to the use of more recent data in our case. 

4.2 @ 14 TeV 

Next, we present a projected study of for the 14 TeV run at 3000 fb“^ of integrated 
luminosity. It should be noted here that the systematic uncertainties used here are for the 
8 TeV run and we have assumed that they will not change significantly for the HL-LHC at 
14 TeV. The inner bands, more clearly noticeable in Fig. 4b, contain only the uncorrelated 
theoretical errors, while the statistical and systematic errors are compounded in the outer 
bands. Clearly, the uncertainty gets reduced, as compared to (@ 7 + 8 TeV), and we 
get an even smaller window around fww and Jbb ~ [~2.76, —2.65] U [—0.06,0.04] TeV“^ 
as shown in Fig. 4. The difference in this case is that the projected band is around the SM 


- 13 - 









-4 -2 0 2 4 


(TeV-2) 

(a) 



-4 -2 0 2 4 

fww/^^ (TeV-') 

(b) 


Figure 3. (a) TZi versus jww/^ (TeV“^) and (b) same plot in magnified scale. Plots (a) and 
(b) are identical for /bb/A^. The red line is the theoretical expectation in presence of HDOs. The 
inner band (dark green) shows the uncorrelated theoretical uncertainty (UThU) and the outer (light 
green) band shows the total surviving uncorrelated uncertainty (UU) (uncorrelated theoretical + 
statistical + systematic) at 7+8 TeV computed using the and [iww (CMS+ATLAS) results. 
The black dotted line is the corresponding central value. The uncertainty bands correspond to 68% 
CL. 


in contrast to what was shown for the 7+8 TeV case, where the ratio of the experimental 
signal strengths was treated as the reference. 

4.3 7^2 @ 14 TeV 

We now show the potential of 7^2 in deriving bounds on some of the operator coefficients 
at 14 TeV. As is evident from Eq. (3.5), this ratio has the capacity to probe Ovp which 
cannot be constrained from TZi. On the other hand, the operator Obb, though amenable 
to probe via TZi, fails to show any marked effect on 7^2 because gets cancelled in 

the ratio as defined by us. Also, Obb does not modify awH but, 7^2 is however sensitive 
to the operator Oww as both (Tvbf and (Twh are sensitive to this. 

By closely following the ATLAS analyses in the context of high luminosity LHC run, 
we have used a trigger cut of 50 GeV on jet pT, instead of using ? 7 -dependent pT cut for 
jets as used in Ref. [68]. The reason is that, a flat cut on the px will most certainly give us 
a less pessimistic number of final state events than that for the p dependent px cuts and 
performs as good as the r?-dependent cut to suppress the background. So, we estimate a 
slightly larger number of events, i.e. we obtain a better efficiency to the cuts for the flat 
px case as compared to what is predicted by ATLAS. For the WH production mode, we 
use a matched sample with WH + 0,1,2 jets with the W decaying leptonically. Finally 
we demand samples with a maximum of one jet in our analysis. In selecting this 0 + 1 jet 
sample, from a matched two jet sample, we encounter another theoretical scale uncertainty 
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Figure 4. (a) TZi versus ^ww!^ (TeV“^) and (b) same plot in magnified scale. Plots (a) and 
(b) are identical for The red line is the theoretical expectation in presence of HDOs. 

The inner band (dark green) shows the uncorrelated theoretical uncertainty (UThU) and the outer 
band (light green) shows total uncorrelated uncertainty (UU) (uncorrelated theoretical + statisti¬ 
cal -I- systematic) at 14 TeV with 3000 fb“^ integrated luminosity. The black dotted line is the 
corresponding central value. The uncertainty bands correspond to 68% CL. 



-4 -2 0 2 4 

(TeV-') 

(b) 


as described in Ref. [74]. We have estimated this uncertainty as follows: 


l^th. 


a{pp WH + > 2 jets) 


X Aa{pp ^ WH + >2iets){pF,fJ‘R), 

rriH 


(4.2) 


where Aa{pp WH + > 2 jets) is the maximum deviation of the exclusive 2-jet cross- 
section computed at pp = pR = mn from the ones computed by varying pp and pR 
between mp[l2 and 2mii. 

In constructing 7^2) we include the modified cut-efficiencies [9, 14] for both the VBF 
and WH channels. Even though we stick to small values of fi where the modification in 
such efficiencies from the SM-values are small, we still incorporate these to make the study 
more rigorous. In computing the statistical uncertainties, we took the relevant numbers 
from the 14 TeV projected study done by ATLAS (see Refs. [68, 69]). Besides, we also 
suggest tagging a single jet for VRF, which reduces the statistical uncertainty by a factor of 
y/2 [75]. The ^/2 factor takes into account the number of events as well as the contamination 
due to ggF as can be seen on Table 1 in Ref. [75]. In Fig. 5, we present 7^2 as a function 
of the fww and fw taken one at a time for an integrated luminosity of £ = 3000 fb“^. 
The outer band (light green) shows the uncertainties due to the statistical, systematic 
compounded with the uncorrelated theoretical part. The central black dashed line shows 
the SM expectation for 7^2- We can see in Fig. 5 that very small values of HDO coefficients 
can be probed by measuring the observable 7^2- For fww-, one can corner the allowed 
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Figure 5. The ratio TZ 2 versus (a) fww !^ (TeV ^), (b) fw! (TeV for the 14 TeV analysis 
with 3000 fb“^. The red line is the theoretical expectation in presence of HDOs. The inner band 
(dark green) shows the uncorrelated theoretical uncertainty due to PDF+a^, QCD-scale and A*^- 
which is defined in Eq. (4.2). The outer band (light green) shows the uncertainties due to the 
statistical, systematic compounded with the uncorrelated theoretical part. The black dotted line is 
the corresponding SM value. The uncertainty bands correspond to 68% CL. 


region to a small window of [—1.96, +1.62] and for fw the range would be [—2.10, +2.50]. 
Predicting the observability of such small values in the parameter coefficients is definitely 
an improvement on existing knowledge. 

4.4 7^3 @ 14 TeV 

The operator Ob appears only in the HZZ and HZ^ couplings. As seen in Fig. 2b, the 
sensitivity oi Ob is too low and hence H ^ ZZ* will not give a proper bound on /^/A^. 
Recent experiment by ATLAS (CMS) puts bounds on the observed signal strength oi H ^ 
Z"f at about 11 (9.5) times the SM expectation at 95% conhdence level [57, 58]. Instead 
of using these weak signal strengths, we perform an analogous projected study of 7^3 at 14 
TeV in the same spirit as at 14 TeV. From Fig. 6 , we find that the projected bounds on 
/b/A^ is [—8.44, —7.17] U [—0.72, +0.56]. The region in between is again inaccessible due to 
poor statistics, as in this region, becomes insignificant, the reasons being similar 

to those mentioned for H —?> 77 . The inner band (dark green) includes the uncorrelated 
theoretical uncertainties due to the partial decay widths of iL —)■ Z 7 and H —)• WW*. The 
outer band (light green), in addition to the theoretical uncertainties, contains the statistical 
and systematic uncertainties. As discussed earlier, a few types of correlated systematic 
uncertainties related to the uncertainty in luminosity, lepton identification and isolation 
etc. will get cancelled in the ratio TZ^. On the other hand, photon identification, isolation 
etc. uncertainties will retain in the analysis. In Table 7, we summarize our obtained region 
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/b/A 2 (TeV-') 


(a) 


(b) 


Figure 6. The ratio TZ^ versus fs !(TeV at 14 TeV with 3000 fb The red line is the 
theoretical expectation in presence of HDOs. The inner band (dark green) shows the uncorrelated 
theoretical uncertainty (UThU) and the outer band (light green) shows the total uncorrelated uncer¬ 
tainty (UU) due to statistical, systematic and the uncorrelated theoretical part. These uncertainty 
bands are for 7?,3 at 14 TeV. The black dotted line is the corresponding SM value. The uncertainty 
bands correspond to 68% CL. 


Observable 

Oww 

Obb 

Ow 

Ob 


[-3.32,-2.91] 

[-3.32,-2.91] 

Not 

Not 

TZi @ 7+8 TeV 

U 

[+0.12,+0.57] 

U 

[+0.12,+0.57] 

bounded 

bounded 


[-2.76,-2.65] 

[-2.76,-2.65] 

Not 

Not 

TZi @ 14 TeV 

U 

[-0.06,+0.04] 

U 

[-0.06,+0.04] 

bounded 

bounded 

772 @ 14 TeV 

[-1.96,+1.62] 

Not 

bounded 

[-2.10,+2.50] 

Not 

bounded 


Not 

Not 

Not 

[-8.44,-7.17] 

773 @ 14 TeV 

used 

used 

used 

U 

[-0.72,+0.56] 


Table 7. We summarize our obtained allowed region of the coefficients of HDOs using the three 
observables. TZz is not used to constrain the operators Oww^Obb and Ow as has been discussed 
in Sec. 3.3. 

of the parameter space allowed using three ratios, TZi, 7^2 and TZ^- We present TZi using 
combined ATLAS-I-CMS data for 7-1-8 TeV run. We also present a projected study for all 
three observables at 14 TeV with an integrated luminosity of 3000 fb“^. The allowed regions 
on fww and fsB shrink at the 14 TeV 3000 fb“^ run as compared to the current data. 
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Using the ratio, 7^2 one can also put bounds on fww and fw- As mentioned earlier, there 
is a ‘special’ region of parameter space where TZi mimics the SM expectation, therefore, 
7^2 can also be used to infer the presence of Oww with ‘special’ values of coefficient fww- 
The operator Ob does not show any appreciable sensitivity in any production of Higgs or 
its decay except in the Therefore, the ratio TZ^ is constructed to constrain fs 

by a significant amount as evident from Table 7. 

5 Summary and conclusions 

We have investigated how well one can constrain dimension -6 gauge-invariant operators 
inducing anomalous HVV interactions. Probing the gauge invariant operators individu¬ 
ally, we feel, are important, since they can point at any new physics above the electroweak 
symmetry breaking scale. While the operators contributing to 77 —)• 77 are subjected to 
the hitherto strongest limits using the (7-1-8) TeV data, the remaining ones are relatively 
loosely constrained, in spite of the bounds coming from precision electroweak observables. 
At any rate, it is necessary to reduce uncertainties as much as possible, since any realisti¬ 
cally conceived new physics is likely to generate such operators with coefficients no greater 
than Ri 0(1) TeV“^. We show that a good opportunity to probe them at this level, and 
improve spectacularly over the existing constraints, arises if event ratios in various channels 
are carefully studied. These include both ratios of events in different final states with the 
same Higgs production channel and those where a Higgs produced by different production 
modes ends up decaying into the same hnal state. While a majority of the theoretical 
uncertainties cancel in the former category, the latter allow us to probe those cases where 
some dimension -6 operators shift the rates in the numerator and the denominator in op¬ 
posite directions. We find that, after a thorough consideration of all uncertainties, all the 
couplings can be pinned down to intervals of width ~ 0(1) TeV“^ on using 3000 fb~^ 
of integrated luminosity at 14 TeV. Even with 300 fb“^, the improvement over existing 
constraints is clearly expected, and the results are more uncertainty-free than in any other 
hitherto applied method. However, we must mention here that this approach should be 
complemented with the study of differential distributions which is not within the scope of 
this paper. 
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